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Abstract 



The atmospheres of Jupiter and Saturn exhibit strong and stable zonal winds. How deep the winds 
penetrate unabated into each planet is unknown. Our investigation favors shallow winds. It consists 
of two parts. 

The first part makes use of an Ohmic constraint; Ohmic dissipation associated with the planet's 
magnetic field cannot exceed the planet's net luminosity. Application to Jupiter (J) and Saturn (S) 
shows that the observed zonal winds cannot penetrate below a depth at which the electrical conduc- 
tivity is about six orders of magnitude smaller than its value at the molecular-metallic transition. 
Measured values of the electrical conductivity of molecular hydrogen yield radii of maximum 
penetration of 0.96Rj and 0.86i?s, with uncertainties of a few percent of R. At these radii, the 
magnetic Reynolds number based on the zonal wind velocity and the scale height of the magnetic 
diffusivity is of order unity. These limits are insensitive to difficulties in modeling turbulent con- 
vection. They permit complete penetration along cylinders of the equatorial jets observed in the 
atmospheres of Jupiter and Saturn. 

The second part investigates how deep the observed zonal winds actually do penetrate. As 
it applies heuristic models of turbulent convection, its conclusions must be regarded as tentative. 
Truncation of the winds in the planet's convective envelope would involve breaking the Taylor- 
Proudman constraint on cylindrical flow. This would require a suitable nonpotential acceleration 
which none of the obvious candidates appears able to provide. Accelerations arising from entropy 
gradients, magnetic stresses, and Reynolds stresses appear to be much too weak. These con- 
siderations suggest that strong zonal winds are confined to shallow, stably stratified layers, with 
equatorial jets being the possible exception. 



1. Introduction 



Jupiter and Saturn, are composed primarily of hydrogen and helium with small additions of 
heavier elements. Their atmospheres e xhibit strong, st able zonal winds composed of multiple jets 



associated with azimuthal cloud bands (|Ingersoll 



1990). Zonal winds peak in the equatorial region 
reaching ~ 100 m s _1 on Jupiter a nd ~ 400 m s" 1 on Sa turn[j] The latitudes of Jupiter's jets have 



not changed for at least 8 years (ISmith & Hunt 



within 10% over 25 years (|Porco et al 



2003) 



1976) and their velocities have been constant 



The depth of the zonal winds is unknown. Both deep and shallow flow models have been 
proposed. Wind speeds measured by the Galileo probe at 7.4°iV on Jupiter rose fr om 90 m s" 1 at 



0.4 bar to 1 80 m s 1 at ~ 5 bar and then remained nearly constant until 22 bar (|Atkinson et al. 



1997 . 



1998). It is important to bear in mind that these measurements only sample the winds in 



the outer 1% of the planet's radius where the electrical conductivity is low. In regions of high 
electrical conductivity, the magnetic field lines are frozen into the fluid. Winds in these regions 
woul d cause changes in th e external magnetic field. By comparing Galileo and Pioneer/Voyager 



data, 



Russell et al 



(|2001 alibi) find that increases of 0.3 deg in the dipole tilt and 1.5% in Jupiter's 



dipole moment may have taken place between 1975 and 2000. The form er could be account ed for 



by meridional flows with speeds in the deep interior of order 0.1 cm s 1 (|Guillot et al 



2004) 



Busse ( 



( Tavloi 



1976 



1983. 



19941) advocates deep flows. He applies the Taylor-Proudman theorem 



1923|) to deduce that zonal flows extend along cylinders oriented parallel the rotation 



axis in the molecular envelope, and then terminate at the outer boundary of the metallic core 
where he assumes that hydrogen undergoes a first order phase transition. But data from shock 
wave experiments show that hydrogen undergoes a continuous transition from a semi-conducting 



Wind speeds on Jupiter are determine d relative to System III coordinates which rotate with the angular speed of 



the planet's magnetic field (IDesslei 



19831) . Only differences among wind speeds on Saturn are known because the 



planet's internal rotaton rate is uncertain. 
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molecular state to a highly conducting metallic state as the pressure increases. This contradicts 
the assumption of a first order phase transition at the core-envelope boundary. 



Recently, a modified deep flow model for Jovian zonal flows has been proposed based on 



simulations of conve ction in a shell with a lower boundary near 0.9Rj (|Aurnou & Heimpel 



Heimpel et al. 



2004; 



2005[). The physical meaning of the lower boundary in the modified 



deep flow 



2004 ). So 



model is obscure. Hydrogen cannot undergo a phase change at that radius (|Guillot et al. 
how might the Taylor-Proudman constraint be violated in order to reduce the zonal flow to a near 
zero value below that boundary? We address related issues in §(5] 



In shallow flow models, the observed high-speed flow is confined to a thin, baroclinic layer 
near the cloud level; the interior flow is much slower. Even if the high velocity flow is confined to a 
shallow layer, its forcing may occur at depth. For example, if the flow were to arise from a process 
that conserved angular momentum per unit volume, pU would be approximately conserved, where 
p is the density and U is the magnitude of the flow velocity. Since the density in the interior is 
several orders of magnitudes larger than that near the surface, the flow velocity could then be much 
greater near the surface. On the other hand, the observed zonal flow might be generated by shallow 
forcing due to the turbulence injected at the cloud level by moist convection, differential latitudinal 
solar heating, latent heat rel ease from condensation of water, or other weather layer processes 



(IVasavada & Showman 



20051) . From the thermal wind equation, a latitudinal temperature gradient 



of about 5 — 10K across a few pressure scale heights below the cloud level would cause substantial 
vertical shear, which makes the flow velocity much greater near the surf ace than deeper down 



dlngersoll & Cuzzi 



1969 



Ingersoll et al. 



19841 : 



Vasavada & Showman 



2005). 



The plan of our paper is as follows. Relevant details of the electrical conductivity of molecular 
hydrogen as measured in shock wave experiments are presented in §[2l Sections [3] and |4] are 
devoted to the calculation of Ohmic dissipation based on the assumption that the zonal wind 
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penetrates the planet along cylinders. In the former, the poloidal magnetic field is determined 
by downward extrapolation of the external field. This procedure is appropriate in regions where 
the magnetic Reynolds based on the convective velocity, R c m <C 1 . In the latter, we examine the 
consequences of assuming that the poloidal magnetic field is parallel to the rotation axis in regions 
where R c m 3> 1. The requirement that the total Ohmic dissipation be bounded from above by 
the planet's net luminosity, C, limits the depth to which the observed zonal winds can penetrate. 
Section [5] asks whether the zonal winds might be truncated within the convective envelope. A 
short summary of our main results is given in §[6] A few technical details are relegated to the 
Appendix. 



2. Electrical Conductivity In Jupiter & Saturn 

Electrical conductivity in the interiors of Jupiter and Saturn is due mainly to hydrogen. Near their 
surfaces it might be significantly enhanced relative to pure hydrogen by the addition of some more 
readily ionized heavier elements. Helium is unimportant due to its high ionization potential. 



Condensed molecular hydrogen is a wide band-gap insulator at room temperature and pressure, 
with a band gap, E g , of about 15 eV, corresponding to the ionization energy of the hydrogen 
molecule. As the pressure increases, this gap is expected to diminish and finally close to zero, 
resulting in an insulator-to-metal transition. In experiments, this transition appears to be gradual. 
As the energy gap closes, hydrogen molecules be gin to dissociate t o monatomic hy drogen 



and electrons start to be delocalized from ions dNellis et al 



1996 



Weir et al. 



1996). The 



insulator-to-metal transition is expected to occur even though the hydrogen molecules have not 
been fully pressure-dissociated. At much higher pressure and temperature, molecular dissociation 
become s complete and it is presumed that pure monatom ic hydrogen forms a metallic Coulomb 



plasma (|Stevenson & Ashcroft 



19741 : 



Hubbard et al 



1997) 
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The conductivity of hydr ogen has been measured in reverb erating shock wave experiments 



from 0.93 Mbar to 1.8 Mbar (Weiretal 



from 0.1 Mbar to 0.2 Mbar (NeJlisetaL 



1996 



Nellis et al 



1999) and in single shock experiments 



1992). In these experiments, hydrogen is in thermal 



equilibrium at pressures and temperatures similar to those in the interiors of giant planets. From 
0.93 to 1.4 Mbar, the measured electrical conductivity of hydrogen increases by four orders of 
magnitude. Between 1.4 Mbar and 1.8 Mbar, the conductivity is constant at 2 x 10 5 S m -1 , similar 
to that of liquid Cs and Rb at 2000 K and two orders of magnitude lower than that of a good metal 
(e.g. Cu) at room temp erature. The cons tant conductivity suggests that the energy gap has been 



thermally smeared out (|Weir et al. 



19961) . Temperature s of shock-compressed liquid hydroge n 



1995 



Holmes et al 



1995) 



have been measured optically in separate experiments (INellis et al 
At the highest obtained pressure of 0.83 Mbar, the measured temperature of 5200 K falls below 
that predicted for pure molecular hydrogen. This is due to the dissociation of molecular hydrogen 
and enables us to estimate the fractional dissociation as a function of pressure. At 1.4 Mbar 
and 3000 K, the dissociation fraction is ~ 5%. Thus metallization of hydr ogen occurs in the 



diatomic molec ular phase and is caused by electrons delocalized from ions (jNellis et al. 



Ashcroft 



1996; 



1968|). Since we are interested in the outer shell of the giant planets, the measurements 



at 0.1 — 0.2 Mbar are the most relevant. In this low-pressure range, the dissociation of hydrogen 
molecules is unimportant. 



The electrical conductivity of a semiconductor can be expressed in the form: 

E 9 ( P ) 



a = (ToO)exp 



2K B T 



(1) 



where a is electrical conductivity, E g (p) is the energy of the density dependent mobility gap, Kb 
is Boltzmann's constant, T is the temperature, and exp {—E 9 /2KbT) expresses the fractional 
occupancy of the curren t carrying states. The conductivity measurement at 0.27 g cm~ 3 and a 



temperature of 4160 K PI (jNellis et al 



(IGuillot 



19921) is close to the interior isentropes of Jupiter and Saturn 



1999J). This measurement determines a = 1.1 x 10 8 S m -1 and E g = 11.7 ± l.leV. The 



2 The corresponding pressure is 0.187 Mbar. 



error bar in the energy gap comes from the experime ntal uncertainties in a and uncertainties in 
calculated post-shock temperatures (jNellis et al.ll 19951) . 



The data suggest that it is sufficiently accurate to assume E g i s a linear function of hyd rogen 



density p H2 : E g = a + bpn 2 , where a and b are two constants (INellis et al. 



1992, 



19961) . The 



measurements and the energy gap (15eV) at room temperature and pressure determine the values 
of a and b. Since the volume mixing ratio of hydrogen in the outer shells of giant planets is 
about 92%, the density in gi ant planet int eriors is 1.18p# 2 . With a , E g (p H2 ), and T(p) from 



the planetary interior model (IGuillot 
calculated. 



1999), the conductivity profiles of the giant planets can be 



Figure (0Q) displays the conductivity distributions as a function of radius in the outer shells of 
Jupiter and Saturn. The conductivity increases with depth, and there is a smooth transition from 
semi-conducting to metallic hydrogen at pressure (1.4 Mbar). This transition takes place at about 
0.84 of Jupiter's radius and 0.63 of Saturn's radius. 



The electrical conductivity is proportional to the total number density of electrical charge carri- 
ers: cr oc n e , which includes a contribution, which we have neglected, from impurities x in addition 
to that from hydrogen: 



n. 



71h 2 ex P 



E n 



n x exp 



E x 



(2) 



2K B T J ^ V 2K B T 

' X x 

where n x and E x express the number density of the electrons and the energy gap due to an 



impurity. Alkali metals are sources of small 



?and gap impuritie s . They may a lso contribute to the 



radiative opacity thus insuring adiabaticity (IGuillot et al. 



2004; 



Guillot 



20051) . The mixing ratio 



of an alkali metal in the interior of a giant planet is presumably similar to that determined from 
its cosmic abundance. With these abundances, a band gap of a few electron volts would lead to a 
conductivity of 10~ 6 ~ 10~ 4 S m" 1 at T ~ 1000 K, significantly above the value due to hydrogen 
in the outer shells of giant planets. 



In magnetohydrodynamics it is conventional to characterize the electrical conductivity a in 
terms of the magnetic diffusivity A = (/ioa)™ 1 , where /i is the magnetic permeability. Figure 
(OQ) shows that the electrical conductivity of hydrogen decreases exponentially outward from the 
metallic conducting region. Therefore, the magnetic diffusivity increases exponentially outward. 
The scale height of magnetic diffusivity 

is shown in figure ©• 



3. Ohmic Dissipation Based On Inward Extrapolation Of The External Magnetic Field 

We approximate the planet 's magnetic fielc 



and Saturn's less than 0.1° (Connernev 



as axisymmetric; Jupiter's dipole tilt is about 10° 
1993|) . Then we evaluate the azimuthal component of the 
magnetic field produced by differential rotation acting on the poloidal components. The maximum 
penetration depth is that of the level above which the associated Ohmic dissipation matches the 
planet's net luminosity, C. 

To proceed, we need to know the poloidal magnetic field above the maximum penetration 
depth. Here we assume that it can be determined by inward extrapolation of the planet's external 
magnetic field. This assumption is appropriate provided the magnetic Reynolds number based on 
the convective velocity field, R c m , remains small down to the maximum penetration depth, which 
our estimates suggest it doesl^ 

Lack of accurate magnetic field measurements at high latitudes close to Jupiter and Saturn 
makes the inward extrapolation of their external magnetic fields somewhat uncertain. Thus we 
cannot exclude the possibility that where R c m 3> 1 the magnetic field might be closely aligned with 
the rotation axis. This possibility is examined in §@J 



An axisymmetric poloidal field is invariant under differential rotation. 
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a. Derivation ofOhmic dissipation 



The time evolution of the magnetic field satisfies 

<9B 



dt 



V x (U x B) - V x [A(r)V x B] 



(4) 



where U and B denote velocity and magnetic field. We work in spherical coordinates and set 
U = U^e^ = r sin 9 fie^. The generation of toroidal field from poloidal field is described by 

aa* „ /an „ 1 on 



dt 



'dQ 
r sin 9 ( -^—B r 
or 



r 89 



X d 



1 d 
1 d 



r dr I 9r ^ ^ 



.sin 



(5) 



r 2 (96* V sin 9 89 

We seek a steady state solution noting that scales proportional to A" 1 , and H x is much 
smaller than the length scale for the meridional variation of U$ and B. Thus we neglect r~ 1 d/d9 
with respect to d/dr, which is equivalent to assuming that j r <C j^O Then, by integrating the 
steady state version of equation ©, we arrive at 



i R pR 

« — sin 9 / dr'r 

J r 



on 

dr' 



B r 



i on 

^~d9 



B, 



With axial symmetry, 



39 



so 



jo(r,9) 



— sin 9 



R 



dr r 



' '2 



J__d_ 

/j r dr 



on 



7, B r 



1 dVt 



'-Be) + 



RX(R) 
r\(r) 



je(R,9). 



(6) 



(7) 



(8) 



dr' r' 89 

Thus jg(r, 9) is determined up to an unknown function of 9, namely je(R, 9), the current density 
in the ionospherel^ There is a simpler and more intuitive way to derive equation ®. Start from 
V x E = and F = E + U x B. Express the former in terms of a line integral around a circuit 
consisting of two small arcs of the same angular width 59, one at r and the other at R, connected 



A toy problem illustrating the effects that H\ <C H p has on j is presented in the Appendix. 

5 A(r) is extremely large, effectively infinite, in the neutral atmosphere but decreases dramatically in the ionosphere. 
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by radial segments of length 5r = R — r. Eliminate the components of E in terms of those of F 
and U x B. Then use j = crF to replace the components of F in terms of those of j. Make the 
assumption \ j r \ <C \ je \ and equation ® is derived. 

Next, we estimate the size of the rhs of equation ®. Let 7\ and T 2 denote the first and second 
terms, respectively. Steady state, axisymmetric currents cannot close within the ionosphere. 
Nor can they penetrate inside the planet due to the extremely low conductivity of its neutral 
atmosphere. However, they can and do flow along field lines into the magneto sphere. These 
currents produce torques which transfer angular momentum from the planet's spin to plasma 
that is drifting outward in the magneto sphere. In so doing, they tend to maintain that plasma 
in approximate corotation with the planet. The torque, and hence jo(R,9), are proportional to 
VLM, where M is the rate at which plasma is expelled from the magnetosphere. Because the 
zonal winds cause Vt to vary by only a few percent, jg(R, 9) is a weak function of 9. Thus, to a 
first approximation, T x and T 2 are uncorrected, so the magnitude of jg(r, 9) is that of the larger 
of these terms. We take a conservative approach and accept the well-determined value of the 
first term as the minimum value for \je(r, 9)\. This is equivalent to treating the ionosphere as an 
equipotential surface in the reference frame rotating with planet's mean angular velocity. 

Ohmic dissipation per unit volume is equal to the square of the current density divided by 
the electric conductivity. Since jo is dominant, we apply equation © to obtain the total Ohmic 
dissipation above radius r; 

For Vt constant on cylinders, the term in round brackets reduces to {dVL/ dw'^B^, where w = 
rsin9. Thus Ohmic dissipation vanishes both for solid body rotation and for a poloidal field 
aligned parallel to the rotation axis. 
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b. Estimation ofOhmic dissipation inside Jupiter and Saturn. 



We apply equation © to evaluate the tota 



zonal flows observed on Jupiter and Saturn (IPorco et al 



Ohmic dissipation ab ove radius r. Atmospheric 



2003 



2005) are taken to be constant on 



cylinders outside a spherical truncation radius and to vanish inside. Since the observed zonal flows 
are not exactly N-S symmetric, we construct N-S symmetric profiles by reflecting the northern 
hemisphere zonal flow about the equatorQ The magnetic fields of Jupiter and Saturn have been 



measured by various spacec rafts and fit 



?y mo dels dominated by a dipole plus smaller quadrupole 



and octupole components (IConnerney 
models in our calculations. 



19931) . We adopt the axisymmetric part of these field 



Total Ohmic dissipation is plotted against truncation radius in Figure © for Jupiter and Saturn. 
It matches the planet's net luminosity at radii of 0.96 Rj and 0.86 R$. The magnetic diffusivity 
at the radius of maximum penetration is 10 7 m 2 s _1 for Jupiter and 3 x 10 6 m 2 s _1 for Saturn. 
By comparison, the magnetic diffusivity is at about 4 m 2 s _1 at the planet's outer metallic cores 
located at 0.84R/ and 0.63i?s, respectively. 

The magnitudes of the induced toroidal magnetic field and the associated poloidal current are 
each inversely proportional to A and thus increase inward. In figure © we display the toroidal 
magnetic field as a function of co-latitude at the maximum penetration depth. It reaches a magni- 
tude of about 15.0 G for Jupiter and about 0.5 G for Saturn. 

The above estimates are based on the downward continuation of the observed poloidal mag- 
netic field. This is a reasonable procedure provided the magnetic Reynolds number based on the 
convective velocity is much smaller than unity above the radius of maximum penetration. 

6 We have verified that using the reflected southern hemisphere zonal flow makes a negligible difference to our 
results. 
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4. Ohmic Dissipation For Poloidal Field Lines Aligned With The Rotation Axis 

Magnetic fields produced in some dynamo simulatio ns show significant alignment of poloidal field 



lines parallel to the rotation axis (|Glatzmaiei 



2005). Such alignment would reduce Ohmic dissi- 



pation associated with differential rotation since this dissipation is proportional to [dQ,/ l dw)B VJ . 
However, alignment can only occur in regions where R c m > 1, so its overall effect on the maximum 
penetration depth of atmospheric zonal winds on Jupiter and Saturn is not obvious. 

To examine the effects of alignment on Ohmic dissipation, we consider a model in which the 
magnetic field is perfectly aligned inside a sphere of radius < R (See figure |5]), 

B m = and B z =B [l- {w/r*f] , (10) 

with p a positive integer. The constant B is set to match the planet's external magnetic dipole. For 
larger p, the internal magnetic field is more concentrated towards the rotation axis. 

The spherical shell r = r* marks the outer boundary of the dynamo region. For r < r*, 
R c m > 1 and the poloidal components of the fluid motions and magnetic field are strongly coupled. 
We assume that the dynamo maintains the aligned poloidal magnetic field against Ohmic decay. 
For r > r*, the poloidal magnetic field is taken to be a potential field. Differential rotation acting 
on the poloidal magnetic field produces a toroidal magnetic field as described in §0 

Outside r„VxB = fl, 

oo 

B r = J2 2n ( R / r ) 2n+1 A°„-l(cOS0)<72n-l 
n=l 

oo 

B e = -^(iVr) 2n+1 P 2 Vi(cos% 2n _!. (11) 

71=1 

In order to match the internal field with the external field, we expand the internal field into 
spherical harmonics at r*, 

oo 

B r = £oX^2n-l^-l(cOS0) 

n=l 

oo 

B e = B J2^2n-lP2n-l(^s9) , (12) 



n=l 
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with 

V2n-1 = ~ ^ I 



(4n- 1 

w 2n -i = 



2 J dxx P P^n-ii 

i( 2n-hL dXX ^ (1 - X2)V2PL - Ax) - <B) 



The radial component of the magnetic field is continuous across r*. Thus 

= {r*/Rf n+1 v 2n ^B /(2n) , (14) 

where B = 2(R/r :t ) 3 g 1 /v 1 in order that g 1 match the planet's external dipole moment. 

Ohmic dissipation comes from three sources: the surface current at r#, the current associated 
with the nonuniform field inside r* (for p > 1), and the current which arises from interaction of 
the vacuum field with the zonal flow outside r*. We treat each of these in turn. 

a. Dissipation due to surface current 

Because Bg is discontinuous across r*, the associated surface current, J s = ABg/ /j, , would give 
rise to infinite Ohmic dissipation. However, the transition between internal and external field 
should be spread across a length scale of order the scale height of the magnetic diffusivity, H\. 
Then Ohmic dissipation from the surface current evaluates to 

P~~FJ / dx{AB e (x)) = —r. rr r- > — + W 2n -l) ■ (15) 

Plots of Ohmic dissipation due to the surface current as a function of r* for different values of p 
are displayed in figure ©. 

b. Dissipation due to current inside r*. 
The azimuthal current density 

= ^ffl, = ?M*M/ 1 .^ (16) 



Ho dw no \r 2 J \ rl / 

produces Ohmic dissipation in a layer of thickness H\ given by 

po Jo \ozu J (4p- 5)(4p- 7)ho 

Figure © displays the internal dissipation as a function of r* for different values p. 
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c. Ohmic dissipation due to action of the zonal wind on the poloidal magnetic field. 

7^ for r > r*, so Ohmic dissipation results from the action of the zonal wind on the poloidal 
magnetic field. We evaluate it with the aid of equation © and plot the results in figure ([8]). 

d. Total Ohmic dissipation 

The total Ohmic dissipation from all three sources is plotted as a function of r* in figure © for 
different values of p. Comparison with figure © reveals that alignment does not significantly in- 
crease the maximum penetration depth of zonal winds on Jupiter and Saturn beyond that calculated 
in §|3]by downward extrapolation of their external magnetic fields. 

The arguments in this section do not refer to the magnetic Reynolds number, R c m , associated 
with the convective velocity. Thus the maximum penetration depth we deduce is independent of 
this quantity. Nevertheless, in reality, the model field we investigate only makes physical sense if 
R c m > 1 at TV 

5. How Deep Do The Zonal Flows Penetrate? 

Together, §[3] arid §@]place an upper limit on the depth to which the zonal winds observed in the at- 
mosphere of a giant planet could penetrate. Here we address the more difficult issue of how deeply 
they actually do penetrate. In particular, do they extend unabated into the convective envelope. It 
proves convenient to work in cylindrical coordinates (w, (ft, z). 

The Navier-Stokes equation which governs the motion of the fluid reads 



where U is the total velocity measured relative to a frame rotating at angular speed Q P (the as- 
sumed uniform angular velocity of the planet's metallic core), p is the density, p is the pressure, 
and $ tot is the gravitational plus centrifugal potential!] We express the steady-state limit of this 

7 We ignore the utterly negligible viscous stress. 



<9U 
~dt 



+ (U ■ V) U + 2Q P x U 




(V x B) x B 



(18) 
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equation as 



— Vp — 



(19) 



where the overbar denotes time average and 



^< VXB)XB -(ITV)U. 

PoP 



(20) 



Smaller terms that involve fluctuations of p, p, and $ tot have been discarded. Individually, the first 
and second terms on the rhs of equation (fl9l) are larger, by far, than the others. Thus 

(21) 



5 - -v* tot 

P 



Taking the curl of equation < fT9l ), we arrive at 

dUfj, d In p 



(9; 



ds 



(Vs x g) • - (V x a) ■ e 9 



(22) 



with g = — V$tot- Thus dU^/dz = for an uniform composition isentrope in the absence of 
Reynolds and Maxwell stresses. This is the Taylor-Proudman state. Next we bound departures 
from this state that each term on the rhs of equation (1221) could produce. 



a. Buoyancy 

1). In the con vective envelope 

Here we are entering uncertain territory. Convection transports the net luminosity in the fluid 
interior, but our understanding of turbulent convection is limited even for nonrotating systems. 
Rotation and especially strong differential rotation add additional complexity. It is not obvious 
which of the terms on the rhs of equation (122)) is dominant for conditions appropriate to a convective 
envelope. 



We are guided by 



Ingersoll & Pollard! (|l982l) who model convection under conditions of strong 



differential rotation which shears convective cell in the azimuthal direction. They argue that under 
these conditions, the magnitude of the component of Vs along g must satisfy 
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dlnp 



ds 



|Vs • g| 



dw 



(23) 
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In other words, the Richardson number based on the rate of shear is of order unity. 

The variation of XJ $ along z depends upon the component of Vs that is orthogonal to both g 
and about which mixing length models are silent. We parameterize this component in terms of 
Vs • g and the angle 5 between Vs and g. Buoyancy drives convection, so it is to be expected that 
8 -C 1. Combining equations (|22)) and (|23T) yields 

2 



HpdUi 






Utj, dz 




ydw J 



tan 5 . 



(24) 



Since g is approximately aligned along the spherical radial direction, |Vs| tan S is essentially the 
magnitude of latitudinal component of Vs. 



At the maximum penetration depth, 



H p dU 4 



for Jupiter, and 



Us dz 



H p dU 4 



0.6 tan 5 



0.3 tan 5 . 



(25) 



(26) 



dz 

for Saturn. As stated above, we expect that tan 5 1. Moreover, its numerical coefficients is 
proportional to H p , which decreases outward. Thus buoyancy is unlikely to effect a significant 
departure from the Taylor-Proudman state in the convective envelope. 



2). In the radiative atmosphere 

There are conflicting views about the strength of the static stability in the radiative atmospheres 
of Jupiter and Saturn and the depth to which stable layers extend. Here we consider how Taylor- 
Proudman columns might be truncated in strongly stable layers in which s changes on the same 
spherical radial scale as p. In such layers 

<9 hip 



in equation (1221) to obtain 



ds 



H p dU 4 



H p \Vs\ 



U 6 dz 



(J 



2VL P U n 
14 



=^ tan 5 . 



(27) 



(28) 



Numerical evaluation yields 



for Jupiter, and 



H p dU 4 



U 6 dz 



H p dU 4 



600 tan 5 



200 tan 5. 



(29) 



(30) 



U4, dz 

for Saturn. Taylor-Proudman columns could be truncated in the radiative atmosphere provided 
surfaces of constant entropy were inclined to those of constant potential in the latitudinal direction 
by more than a fraction of a degree. The resulting velocity field would constitute a strong thermal 
wind. 



b. Magnetic stresses 

Deviations from the Taylor-Proudman state caused by Maxwell stresses follow from 

e,* ■ V x 



HpdU, 



Us dz 



2fi Q P U7 9 



(V x B) x B 



P 



H P B 2 



2fi Q P \U^\p£ 2 



(3D 



where £ is the typical scale over which B varies. The order of magnitude estimate assumes i < H p . 
Next we bound B 2 /£ 2 by applying the Ohmic constraint. The latter reduces to 



4ir\R 2 H p B 2 < 

Eliminating B 2 / £ 2 between equations (1311) and (1321) . we arrive at 

_C 

8n\tt P \U^\pR 2 ' 



(32) 



H p 8U q 



U A dz 



< 



Numerical evaluation yields 



H p dU 4 



<10" 



(33) 



(34) 



U4 dz 

at the maximum penetration depths in Jupiter and Saturn. Moreover, this ratio decreases sharply 
outward. Clearly, magnetic stresses are incapable of truncating the observed zonal flows. 
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c. Reynolds stresses 



Departures from the Taylor-Proudman state caused by Reynolds stresses obey 



H p dU 4 



U A dz 



2VL P \U < , 



e 6 - Vx ((U- V)U 



(35) 



We treat separately contributions from the mean axisymmetric flow, U, and from fluctuations about 
it, u = U - U. 

The sole contribution involving U& can be absorbed by adding U^/w to Q P on the lhs of 
equation (fT9b and in what follows. We do not consider it further. The poloidal part of U, denoted 
by U p , describes meridional circulation. It yields 



H p dU s 



U A dz 



u 2 

p 



(36) 



2{lp\U ( j > \Hp 

where the variation scale of U p is set to H p . With parameters appropriate to Jupiter and Saturn, 
the rhs of equation (l36l) is <C 1 for \U P \ <C \Us\- 
Velocity fluctuations contribute 



H p dU 4 



U d dz 



(37) 



2tt P \U (f) \£ p £ (l> 

where £ p and 1$ are the meridional and azimuthal scales of the convective cells. In deriving equa- 
tion (T371) . we note that that azimuthal stretching of eddies by strong differenti al rotation results in 
'■ 'a > £ p and mass conservation implies u p £s ~ u^£ p . To bound u 2 ,, we follow llngersoll & Pollard 



(|1982|) and adopt \dU 'dm\(£ p / '£</,) as the convective rate. In each scale height, turbulent mechan- 



ical energy is dissipated as heat at a rate ~ C. Hence 



A<KR 2 H-- 2dU ^ p 



C. 



(38) 



Together, equations (1371) and (1381) yield 



H p dU q 



Us dz 



C 



< 



(39) 



where A9 denotes the typical latitudinal width of the zonal jets. In accord with 



Ingersoll & Pollard 



(|1982|) . we set \dUs/dw\ ~ \Us\/£ P in arriving at the final form of equation (l39l) . Numerical 
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evaluation with parameters appropriate to the tops of the convection zones in Jupiter and Saturn 

gives 

< 10" 5 . (40) 



HpdU^ 



U4, dz 

d. Maximum width of an equatorial jet 

A sufficiently narrow equatorial jet could maintain constant velocity on cylinders throughout the 
planet. As an example, we consider the specific velocity profile 

— (tx (6-60) 

U^ = U^sm - / /o °M if 6<n-6 ; (41) 



2(tt/2-0 o ), 
and, 

U (j> = if 6 < 6 and 9 > vr - 6 ; (42) 

so the jet has equatorial velocity U^o and angular half-width n/2 — 9q. For Jupiter and Saturn, 
U^o is approximately 140 m s^ 1 and 400 m s -1 , respectively. Figure (flOl) displays the calculated 
Ohmic dissipation rate as a function of the jet half- width. The maximum half- width is about 18.5° 
for Jupiter, and 35.5° for Saturn. The maximum half widths are related to the fractional radii of 
maximum penetration, r mp /R, calculated in section[3]by 

cosflo*^. (43) 



6. Discussion. 

The condition that the total Ohmic dissipation not exceed the planet's net luminosity sets an upper 
bound on the depth to which the zonal flows observed in the atmospheres of Jupiter and Saturn 
could penetrate. At these depths, the magnetic Reynolds number, based on the observed zonal 
winds and the scale height of the magnetic diffusivity, is of order unity. 

We consider it unlikely that the observed flows extend to the depth of maximum penetration 
because that would require the Taylor-Proudman constraint to be violated in the convective en- 
velope. We have been unable to identify any plausible mechanism that could do this. Lack of a 
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ysis of the robustne s s of th e Taylor 



Ingersoll & Pollard ( 



1982h to apply 



rigorous model for convection is a weakness. We base our ana 
Proudman constraint on the mixing length model designed by 
under conditions of strong differential rotation. However, we have confirmed that the conclusions 
we draw from it also follow from other convection models that consider only solid body rotation, 
or even those that neglect the effects of rotation entirely. 

The boundaries of the cylindrical extensions of the equatorial jets essentially coincide with the 
maximum penetration depths. Thus these jets could maintain constant velocities along cylinders 
through the planets. Winds measured at 7.4°iV by the Galileo probe as it descended into Jupiter 
are consistent with this possibility. 

Acknowledgments We are grateful for support from the NASA Planetary Geology and Geo- 
physics program. 

7. Appendix 

We calculate the current distribution arising from a vacuum axial dipole field and a simple zonal 
flow. The angular velocity, f2, is taken to be of the form 

n = n P + — — S*z* -. (44) 

1 + exp [c mp {r mp - r)\ 

Here e mp is the .rnncaiion fac.orand r mt is the truncation radins. «(,) « P for r < r mv and 
Q(r) ->• Sip + Vt out for r > r mp J3 

The current density arising from interaction of this flow with a vacuum dipole magnetic field 
satisfies: 

V x 1 = V x {[(n - Sl P ) e z x r] x B} = -^^^ , (45) 
a 3r z ar ad 

where M is the dipole moment. 

Because V • j = 0, j can be derived from a vector potential A such that 

j = V x , (46) 



8 We plot fl — Hp as a function of scaled radius in figure ( TfTT i. 
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where 



By direct computation 



a 



6ra 



A 



-rf(r) dP<i 
~~6 d9 



6/ 



d In a d 

dr dr 



(r 2 f) 



de 



So 



— r 



dn 

dr ~ AaM 



dr 2 dr dr 



> 2 f) 



(47) 



(48) 

(49) 
0; Near 



At the planet's surface, the current density along the radial direction goes to zero: f(r) 
its center, j oc r: df(r)/dr — f(r)/r — > 0. We solve equation (|49l to obtain f(r). 
Since, 

j • V (Ar sin e) = , (50) 

the current density follows the contour lines of Ar sin e and its magnitude satisfies | j | = 
I V ( Ar sin e) \ jr sin 6. Figure (fT2l) displays the streamlines of current flow. 
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Current stream lines arising from a toy model. A vacuum dipole field interacts 
with a simple cylindrical flow consisting of a nearly uniformly rotating core and a 
more rapidly rotating envelope. See text and figure (fTTT) for details. In the outer 
envelope, the streamlines are close to lines of constant r. Their spacing is inversely 
proportional to jr sin 9. Ohmic dissipation per unit volume, oc aU^B^, is maximal 
near the dash line which marks to core-envelope boundary. Our model only makes 
sense in regions where the magnetic Reynolds number based on the convective 
velocity is small 
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Figure 1: Electrical conductivity and magnetic diffusivity distributions inside giant planets: (a) 
Jupiter; (b) Saturn. Values of conductivity and magnetic diffusivity are plotted in the left and right 
panels, respectively. Solid lines depict mean value; dashed lines bound the range of uncertainties. 
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(a) Jupiter 



10 6 



(b) Saturn 




Scaled radius Scaled radius 



Figure 2: Scale height of magnetic diffusivity as a function of scaled radius: (a) Jupiter; (b) Saturn. 
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(a) Jupiter 



(b) Saturn 




Figure 3: Observed zonal flow is taken to penetrate along cylinders until it is truncated at r mp . 
Solid curves display nominal values for the total Ohmic dissipation, P, as a function of r mp /R. 
Dashed curves bound its uncertainty. Horizontal solid line marks the planet's net luminosity, C At 

r P = L 
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Figure 4: Toroidal magnetic field at the maximum penetration depth as a function of colatitude: 
(a) Jupiter, (b) Saturn. 
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Figure 5: Poloidal magnetic field lines for p 
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ure 6: Ohmic dissipation due to surface current as a function of r*/R. (a) Jupiter; (b) Saturn. 
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Figure 7: Ohmic dissipation due to internal current as a function of r*/R. (a) Jupiter; (b) Saturn. 
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(a) Jupiter 



(b) Saturn 




Figure 8: Ohmic dissipation due to interaction of £> ro with dVt/dw as a function of r*/R. (a) 
Jupiter; (b) Saturn. 
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Figure 9: Total Ohmic dissipation as a function of r\jR. (a) Jupiter; (b) Saturn. 
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(a) Jupiter 



(b) Saturn 




5 10 15 20 25 30 20 30 40 50 

Angular half-width Angular half-width 



Figure 10: Total Ohmic dissipation rate verses jet's angular half- width. Horizontal solid lines mark 
the planet's net luminosity. Dashed lines indicate the range of uncertainties. Central value of the 
maximum half-width is 18.5° for Jupiter and 35.5° for Saturn. 
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Figure 11: Angular velocity corresponding to current streamlines shown in figure (TT2|) . Parameters 
in equation (|44)l are r mp / R = 0.9, c mp = 50. Surface value of f2p + Vl out is scaled to be near unity. 
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Figure 12: Current stream lines arising from a toy model. A vacuum dipole field interacts with a 
simple cylindrical flow consisting of a nearly uniformly rotating core and a more rapidly rotating 
envelope. See text and figure (fTTj) for details. In the outer envelope, the streamlines are close to 
lines of constant r. Their spacing is inversely proportional to jr sin 9. Ohmic dissipation per unit 
volume, oc aUlBp, is maximal near the dash line which marks to core-envelope boundary. Our 
model only makes sense in regions where the magnetic Reynolds number based on the convective 
velocity is small. 
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